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In this paper, we discuss the c axis optical conductivity 
Re [<7 c (oj)] in the high T c superconductors, in the supercon- 
ducting state. The basic premise of this work is that electrons 
travelling along the c axis between adjacent Cu02 layers must 
pass through several intervening layers. In earlier work we 
found that, for weak inter-layer coupling, it is preferable for 
electrons to travel along the c axis by making a series of in- 
terband transitions rather than to stay within a single (and 
very narrow) band. Moreover, we found that many of the 
properties of the normal state optical conductivity, including 
the pseudogap could be explained by interband transitions. 
In this work we examine the effect of superconductivity on 
the interband conductivity. We find that, while the onset of 
superconductivity is clearly evident in the spectrum, there is 
no clear signature of the symmetry of the superconducting 
order parameter. 
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I. INTRODUCTION 

Recently, there has been a lot of attention paid to the 
c axis optical conductivity in the high T c cuprate su- 
perconductors. In particular, measurements of the opti- 
cal conductivity in YB^CusO^ (YBCOz) have revealed 
that the c axis transport is very different in character 
from the electronicJjjansport in the a and b directions 
(within the layers). En3 There is a great deal of specula- 
tion as to the source of this difference. At one end of the 
spectrum of thought there is the claim that the unusual c 
axis transport is evidence for some ncoi-Fermi liquid like 
ground state within the CuC>2 layers.Q At the other end 
of the spectrum, the claim is that the electronic ground 
state is metallic, but that there is some unconventional 
tunneling mechanism between the layers.u~u The com- 
mon feature to all of these models, however, is that they 
consider only a single copper oxide layer in each unit cell. 

In the case of YBCO^, of course, the situation is not so 
simple. There are several layers between adjacent CuC>2 
layers, one of which (the CuO chap, layer) is known to 
be conducting. In a previous paperEHl we asked the ques- 
tion "What is the effect of these intermediate layers, as- 
suming a simple metallic model?" We considered a sim- 
ple two-layer model in which each unit cell contained a 
CuC>2 plane and a CuO chain, and we calculated the op- 
tical conductivity in the normal state. What we found 
was that, while making almost no assumptions about the 
band structure, we could explain many of the features 
seen in experiment — including the pseudogap seen in op- 



tical experiments. We claim that this work throws into 
doubt any attempt to interpret the c axis conductivity 
that does not take into account the multilayered struc- 
ture of the high T c materials. 

In this work, we examine the c axis optical conductivity 
for a simple model of YBCOj; in the superconducting 
state. This model is the same as the one described above: 
each unit cell contains a plane layer and a chain layer. 
These layers are evenly spaced and connected by coherent 
single-electron hopping. The amplitude for the hopping 
is parameterized by t±. The current model differs from 
our earlier one in two important ways. First, the sample 
is taken to be superconducting. Based on experimental 
observations, we include apSiiperconducting gap in both 
the plane and chain layers.EU Second, the scattering rate 
1 /r is set to zero. In our study of the optical conductivity 
in the normal state, the linear temperature dependence of 
the scattering rate was important at high temperatures. 
At low temperatures the scattering rate is small and can 
be ignored. 

One of the main conclusions of our earlier work is 
that, for small t±, interband processes play a domi- 
nant role in the c axis conductivity. In isolation, the 
plane and chain layers have dispersions £i(k x ,k y ) and 
£,2(k x ,k y ) respectively. When the layers are coupled by 
t±, they hybridize and form two bands e+(k x , k y ,k z ) and 
e-(k x ,k y ,k z ). Electrons which travel along the c axis 
in the presence of an external field may do so by cither 
staying within the bands e±, or by making transitions 
between them. If the layer dispersions £i and £2 are non- 
degenerate, then the bands e± will differ from £1 and 
£> by ~ (2t ± cos(k z d/2)) 2 /(£_ 1 - £>), where d/2 is the 
interlayer spacing (this follows from second order pertur- 
bation theory). The Fermi velocities v z ± — h~ de±/dk z 
therefore scale as t\. On the other hand, the matrix 
element for an interband transition is ~ 2tj_ cos(k z d/2). 
For small t±, the c axis transport occurs preferentially 
through interband transitions. 

While the intraband conductivity has the well-known 
Drude form, the interband conductivity does not. Much 
of this paper will be devoted to understanding how the 
interband conductivity in the plane-chain model differs 
from the usual picture of conductivity in the supercon- 
ducting state. 

We are aware of one other model which attempts to ac- 
count for the presence oiintermediate layers. The model, 
proposed by Abrikosov,E3 examines the effect of resonant 
tunneling through an impurity layer which lies between 
Cu02 layers. In the case of YBCOa;, the impurities are 
the oxygen vacancies in the CuO chain. With his model, 
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he is able to provide a reasonable explanation of the c 
axis d.c. resistivity. 

This paper is organized as follows. In Sec. 0, an equa- 
tion for Re[a c (io)] (the optical conductivity along the c 
axis) is derived. The equation is integrated numerically 
and the results are discussed in Sec. 
sion is contained in Sec. 



IV 



III. A brief conclu- 



II. DERIVATIONS 



The mean field Hamiltonian for our model isQ 

# = E C kMk)C k (1) 



where = 



c ikT' c i-kl>4kT> c 2-kl] and c| k(T creates an 
electron of spin a and 3-dimensional momentum k in the 
plane (i = 1) or chain (i = 2) sublattices. The Hamilto- 
nian matrix is 



h(k) = 



6 A k t(k) 

Ak -6 -i(k) 

i(k) 6 Ak 

-t(k) A k -£ 2 



(2) 



where A k is the mean field superconducting order pa- 
rameter, £i and £2 are the plane and chain dispersions 
respectively, and t(k) couples the plane and chain layers. 
The dispersions £1 and £2 are 

£1 = — 2(7 1 [cos(k x a) + cos(fc y a) — 2B cos(k x a) cos(fcj,a)] — fix , 

(3) 




X ky 



and 



6 = -2cr 2 cos(k y a) - fi 2 (4) 

where a±, a%, /Ui, /12, and B are adjustable parameters, 
and a is the unit cell size in ab directions. For this 
work, the parameters are fixed at {<J\, 02, Mi, M2} = 
{70,100,-65,-175} mcV and B = 0.45. This is done 
by fitting the magnitudes of the penetration depth at 
T = (which effectively measures the magnitudes, of the 
Fermi velocities v x , v y and v z ) to experiment JL3 while 
maintaining a Fermi surface that loaksjqualitatively like 
that of band structure calculations oll3 

The interlayer coupling is t(k) = — 2t± cos(k z d/2), 
which follows from a tight-binding model of c-axis cou- 
pling. The unit cell size is d along the c-axis. Again, 
by fitting A C (T = 0) to experiment, we have determined 
that tj_ - 20 meV for YBC0 6 . 93 and t± - 1 meV for 
YBC0 6 . 7 . 

The Fermi surface is shown in Fig. |l|(a). The dashed 
lines give the Fermi surface of the isolated chain and 
plane layers. These are also the Fermi surfaces of e± 
at k z = n/d, since t(Tr/d) = 0. The solid curves are the 
Fermi surfaces at k z = 0, where t(k z ) is a maximum. For 



FIG. 1. a) Fermi surface for the plane-chain model. The 
dashed lines are the Fermi surfaces of the isolated plane and 
chain layers, the solid lines are the Fermi surfaces at k z — 0. 
For other values of k z , the Fermi surfaces lie between the solid 
and dashed curves, b) Schematic of the energy bands along 
k x — it /a, k z — n/d. The energy bands are shown for the su- 
perconducting state (solid curves). The normal state energies 
e_ and e+ vanish at the points X and Y respectively. Exter- 
nal fields excite two types of interband transition — pair cre- 
ation and quasiparticle transitions — which are shown. In pair 
creation, the final state has one quasiparticle in each of the 
bands. This process is gapped since the energy required for 
the process Tiuj = E+ (k) + E- (k) has a nonzero minimum. In 
quasiparticle transitions, thermally excited quasiparticles can 
then make interband transitions. This process is not gapped 
since there will generally be values of k for which the excita- 
tion energy hui = E+ (k) — E- (k) vanishes. 
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intermediate values of k z , the Fermi surface lies between 
these two curves. 

The superconducting gap Ak is chosen to be the same 
in both the plane and chain layers. It can have either a 
rf-wave symmetry 



A(T)[cos(fc2;a) — cos(fcj,a)], 



or an s-wave symmetry 



A(T). 



In this work we do not propose a microscopic origin for 
the pairing interaction but rather assume that A(T) can 
be described phcnomenologically by 



A(T) 
A(0) 



tanh uw 



T 



with 2A(0)//c B T c = 3.5 and T c = 100.K. 

We have found, in previous work,II9 that Re [er c (u;)] in 
the normal state depends on the choice of parameters in 
two ways: a c (u>) depends qualitatively on whether or not 
the chain and plane Fermi surfaces cross, and on the mag- 
nitude of ij_ . We found that when the Fermi surfaces do 
not cross, there is a direct gap which appears as a pseudo- 
gap of energy Tlujq in a c (u>). The pseudogap disappears 
at high T because of the smearing of the quasiparticle 
energy by the large, temperature dependent, scattering 
rate. The temperature T* at which the pseudogap begins 
to be resolved is therefore determined by 

1/t(T*) ~ cj , 

and there is no direct connection between the onset of 
superconductivity at T c and T* . As we mentioned before, 
the scattering is ignored in this work so the pseudogap is 
unchanged by temperature. We also found in our earlier 
work that for small values of t± the Drude-like intraband 
contribution to the conductivity is hidden by the large 
interband conductivity. 

In this work we will examine the c-axis optical response 
in the superconducting state. We will assume that the. 
sample is in the clean limit. There is good evidencet£l 
that the large, temperature dependent scattering rate 
does, in fact, drop dramatically below T c , although we 
have made this assumption primarily to keep the model 
simple. In this case, the Drude part of a c (uj) collapses 
to a 5-funtion at oj = 0, and the only response at finite 
frequency comes from interband transitions. Since the 
main effect of changing t± is to change the magnitude 
of the interband conductivity and not its form, we can 
arbitrarily fix t± = 20 mcV. 

One of the most interesting features of our model 
Hamiltonian is that the gap Ak is the same in both the 
chain and plane layers. We have discussed the experi- 
mental evidence for this elsewhereO It is a somewhat 
surprising property since the chain layer does not have 
the tetragonal symmetry needed to generate a <i-wave 



order parameter. This strongly suggests that the pair- 
ing interaction must originate in the (tetragonal) CuC>2 
plane layer, but then the difficulty lies in understanding 
the large magnitude of the gap on th&,chain layer. This 

s.tirlHl by Xiang and 



has been jiiscussed elsewhere by us,E 
WheatleyEj and by O'Donovan and Carbotte.Ej 

The eigenvalues of the Hamiltonian matrix /i(k) give 
the band energies E\ = E + , E 2 = —E + , E 3 — £L, 
E± — —E_, where 



£-4 
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and e± are the normal state band energies 

6 + 6 



± 



f 2 . 



(5) 



(6) 



The optical conductivity for a layered system is0 
e 2 h 



Re[cr^(w)] 



E 



dx Tr 



x A(k;ac + ^a;)7 I/ (k,k) 



A(k;x)%(k,k) 
f{x)- f(x + fiio) 



(7) 



where ^4(k; u>) is the spectral function. It is a 4 x 4 matrix 
whose diagonal elements describe the spectral weight in 
the 4 superconducting bands. The electromagnetic ver- 
tex functions 7 for c-axis transport areJlJ 



where 



v z+ aT z 

v z+ (3T Z 

aT z (3T Z v z _ 

-pT z aT z 



-f3T z 
aT z 




A k 



(8) 



E+E. 



1 - 



A£ 



E+E_ 



are the coherence factors and 
1 dt(k z ) 



T' 



h dk z 



6-6 



is the matrix element for the transition. Furthermore v z ± 
are the Fermi velocities 

_ 1 de± 

z± h dk z 

Notice that if the bands £1 and £2 are nondegenerate, 
v z ± cx tj_. In contrast, the off diagonal matrix elements 
in 7 are proportional to Since the diagonal matrix 
elements are for the intraband (or Drude) conductivity, 
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FIG. 2. The conductivity Re [<t c (w)] is shown for a d-wave 
gap and for a range of temperatures. The calculation is in 
the clean limit, so that the conductivity is entirely due to in- 
terband transitions. Notice that even though there is gap in 
the normal state spectrum (at T = T c — 100 K), the super- 
conducting spectrum is gapless at finite temperature. 



we can conclude that, for small t±, the Drude peak will be 
small relative to the non-Drude interband conductivity. 
Further, since a c (uj) depends on -y 2 , the Drude part of Re 
[<7 c (u;)] is proportional to tj_, while the interband part is 
proportional to t\. 

In the clean limit, the spectral function is An(\s.;u>) = 
2it6(uj — Ei), the intraband contribution to a c (u> > 0) 
vanishes and Eq. (0) reduces to 



Re[a c (uj > 0)] 
2f(E 



2ne 2 h 



x < a 



- f(E+) 



E 1 ? 



0- 



E+ -E. 
1 -/(£_)- 



S(Hlo - \Ej, 



E_ 



-5{huj - E^ 



E-\) 
-E-) 



(9) 



This is our basic result. The first and second terms in the 
integrand in (^J) represent interband transitions of ther- 
mally excited quasiparticles and pair creation of quasi- 
particles respectively. We will discuss these two processes 
in more detail in the next section. 



III. RESULTS AND DISCUSSION 

Equation (^) is integrated numerically and the results 
are shown in Figs. || and || for s and ci-wave gaps respec- 
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FIG. 3. The conductivity Re [er c (w)] is shown, as in Fig. 
but for an s wave gap. 



tively. The conductivity is calculated for several temper- 
atures between T = K and T = T c = 100 K. At T = T c 
the system is in the normal state and 



k 



/(C-) "/(£+) 



5(hu) — e + + e_). 



(10) 



Equation (|l^) has the form of a joint density of states 
since the integrand is proportional to 8{huj — e+ + e_). 
The thermal factors ensure that the interband transitions 
are between filled and empty states. At low T, transitions 
occur between states for which e_(k) < and e+(k) > 0. 

In the high T c superconductors, the large frequency 
range over which a c (uj) extends is often ascribed to ei- 
ther a large scattering rate for interplane transitions, or 
to the non-Fermi liquid like nature of the ground state. 
Here, the frequency range is of the order of the bandwidth 
because the energy difference e + (k) — e_ (k) extends over 
a wide range of energies. In Fig. [|(a) the energy differ- 
ence e+(k) — £-(k) is ~ 20 meV at the point X, where 
the Fermi surfaces are close together, and is ~ 375 meV 
at the point Z. 

For our particular choice of parameters, the normal 
state conductivity has a finite band gap, which we iden- 
tify with the pseudogap seen in optical conductivity ex- 
periments. The minimum value of hujQ — e+(k) — e_(k) 
for which a c (uio) ^ is at the point X. Clearly the value 
of loq depends on the distance between the two pieces 
of Fermi surface. If the Fermi surfaces cross, then loq 
vanishes and arbitrarily low energy excitations are pos- 
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sible. In this case there is no pseudogap, and we claim 
that this describes YBCO in the the optimally doped and 
overdoped cases. 

As the temperature is lowered and the material be- 
comes superconducting, two types of interband process 
become possible. As in the normal state, quasiparti- 
cles which occupy one of the bands may make transi- 
tions into the other band. This is described by the first 
term in (Q) and the energy required for the transition 
is fiui = |_E+(k) — _E_(k)|. Below T c , however, most of 
the electrons are in the superconducting condensate, and 
the largest contribution to the interband conductivity 
comes from the creation of quasiparticle pairs, in which 
a Cooper pair is broken and one quasiparticle goes into 
each of the two bands. This is described by the second 
term in (^) , and the energy required for the transition is 
hu = E + (k) + £L(k). In Fig. |, a c {u) is shown as the 
sum of the quasiparitcle and pair creation terms. 

The two types of interband transition are illustrated 
m Fig. 0(b), and we can make a few comments about 
Eq. (Q) just by inspection of the figure. The first is that 
the quasiparticle term is not gapped since excitations of 
arbitrarily low energy are available near k-points where 
E + (k) — _E_(k). The probability for such transitions to 
occur, however, is strongly supressed by the probability 
that the lower energy band is initially occupied. At T = 
0, the thermal factor f(E + ) — /(£"_), ensures that the 
quasiparticle term vanishes. 

A second point we can make is that the pair cre- 
ation term is gapped since the energy of pair creation 
hui = £'+(k) + -E-(k) has, in general, a nonzero mini- 
mum value. The lowest energy pairs which can be created 
can obviously be found by minimizing y/ e + (k) 2 + A k + 
y/e-(k) 2 + Aj[. For the case where the Fermi surfaces 
cross, so that there is a line of k- values along which 
e+(k) = e_(k) = 0, the minimum energy is near Tilo = 
2Ak- Importantly, this shows that the pair creation term 
is gapped for a d-wave superconductor unless the Fermi 
surfaces happen to cross at Ak = 0. 

Perhaps the most striking feature of Figs. and |^ is 
that the symmetry of the gap does not reveal itself in 
any obvious fashion. It is customary — within a Drudc 
model — to associate a gapped frequency dependence with 
an s-wave order parameter, and a gapless frequency de- 
pendence with a d-wave order parameter. It is clear from 
the above discussion, however, that this cannot be done 
here. To summarize our discussion simply, the interband 
contribution to the conductivity — which is dominant at 
low temperatures — probes the structure of the energy 
sum E + + E- and the energy difference | E + — E- \ , and 
not of the gap energy Ak . 

This is in contrast, for example, with the single layer 
model in wduch the c axis transport is through diffusive 
scattering. u u In this model the conductivity probes the 
density of states of a single layer 
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FIG. 4. The interband conductivity is shown at 75 K (solid 
curve) for a d-wave gap. There are two contributions to 
the conductivity: a quasiparticle term results from interband 
transitions of thermally excited quasiparticles (dashed curve) 
and a pair creation term (dot-dashed curve). The quasipar- 
ticle contribution is gapless and vanishes at T = 0. The pair 
creation term is gapped. 



IV. CONCLUSIONS 

For layered superconductors, in which there are more 
than one type of layer, the c axis conductivity is domi- 
nated by interband effects when the interlayer coupling is 
weak. In other words, it is easier for an electron to travel 
along the c axis by making a series of interband transi- 
tions than by staying within a band. The intraband (or 
Drude) and interband conductivities vary as t\_ and t\ 
respectively. 

There are two types of interband transition. The first, 
interband transitions of thermally excited quasiparticles, 
has a gapless spectrum and probes the joint density of 
states of the dispersion E + (k) — E_(k). The second, the 
excitation of quasiparticle pairs, has a gapped spectrum 
and probes the joint density of states of the dispersion 
E + (k) + ^.(k). Since neither of these processes probes 
the structure of a single band, there is no clear signature 
of the gap symmetry in Re [cr c (w)]. 
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